Target of Opportunity Observations of Gravitational Wave Events with LSST

The discovery of the electromagnetic counterparts to the binary neutron star merger GW170817 has opened the era of GW+EM multi-messenger astronomy. Exploiting this breakthrough requires increasing samples to explore the diversity of kilonova behaviour and provide more stringent constraints on the Hubble constant, and tests of fundamental physics. LSST can play a key role in this field in the 2020s, when the gravitational wave detector network is expected to detect higher rates of merger events involving neutron stars (∼10s per year) out to distances of several hundred Mpc. Here we propose comprehensive target-of-opportunity (ToOs) strategies for follow-up of gravitational-wave sources that will make LSST the premiere machine for discovery and early characterization for neutron star mergers and other gravitational-wave sources

A Standard Siren Measurement of the Hubble Constant from GW170817 without the Electromagnetic Counterpart

We perform a statistical standard siren analysis of GW170817. Our analysis does not utilize knowledge of NGC 4993 as the unique host galaxy of the optical counterpart to GW170817. Instead, we consider each galaxy within the GW170817 localization region as a potential host; combining the redshifts from all of the galaxies with the distance estimate from GW170817 provides an estimate of the Hubble constant, H_0. Considering all galaxies brighter than 0.626L*_B as equally likely to host a binary neutron star merger, we find H_0 – 77^(+37)_(-18) km s^(−1) Mpc^(−1) (maximum a posteriori and 68.3% highest density posterior interval; assuming a flat H_0 prior in the range [10, 220] km s^(−1) Mpc^(−1)). We explore the dependence of our results on the thresholds by which galaxies are included in our sample, and we show that weighting the host galaxies by stellar mass or star formation rate provides entirely consistent results with potentially tighter constraints. By applying the method to simulated gravitational-wave events and a realistic galaxy catalog we show that, because of the small localization volume, this statistical standard siren analysis of GW170817 provides an unusually informative (top 10%) constraint. Under optimistic assumptions for galaxy completeness and redshift uncertainty, we find that dark binary neutron star measurements of H_0 will converge as 40%/√(N), where N is the number of sources. While these statistical estimates are inferior to the value from the counterpart standard siren measurement utilizing NGC 4993 as the unique host, H_0 = 76^(+19)_(-13) km s^(−1) Mpc^(−1) (determined from the same publicly available data), our analysis is a proof-of-principle demonstration of the statistical approach first proposed by Bernard Schutz over 30 yr ago

Low Frequency Tilt Seismology with a Precision Ground Rotation Sensor

We describe measurements of the rotational component of teleseismic surface waves using an inertial high-precision ground-rotation-sensor installed at the LIGO Hanford Observatory (LHO). The sensor has a noise floor of 0.4 nrad$/ \sqrt{\rm Hz}$ at 50 mHz and a translational coupling of less than 1 $\mu$rad/m enabling translation-free measurement of small rotations. We present observations of the rotational motion from Rayleigh waves of six teleseismic events from varied locations and with magnitudes ranging from M6.7 to M7.9. These events were used to estimate phase dispersion curves which shows agreement with a similar analysis done with an array of three STS-2 seismometers also located at LHO

Classifying the unknown: discovering novel gravitational-wave detector glitches using similarity learning

The observation of gravitational waves from compact binary coalescences by LIGO and Virgo has begun a new era in astronomy. A critical challenge in making detections is determining whether loud transient features in the data are caused by gravitational waves or by instrumental or environmental sources. The citizen-science project \emph{Gravity Spy} has been demonstrated as an efficient infrastructure for classifying known types of noise transients (glitches) through a combination of data analysis performed by both citizen volunteers and machine learning. We present the next iteration of this project, using similarity indices to empower citizen scientists to create large data sets of unknown transients, which can then be used to facilitate supervised machine-learning characterization. This new evolution aims to alleviate a persistent challenge that plagues both citizen-science and instrumental detector work: the ability to build large samples of relatively rare events. Using two families of transient noise that appeared unexpectedly during LIGO's second observing run (O2), we demonstrate the impact that the similarity indices could have had on finding these new glitch types in the Gravity Spy program

Modeling Multi-Wavelength Stellar Astrometry. I. SIM Lite Observations of Interacting Binaries

Interacting binaries consist of a secondary star which fills or is very close to filling its Roche lobe, resulting in accretion onto the primary star, which is often, but not always, a compact object. In many cases, the primary star, secondary star, and the accretion disk can all be significant sources of luminosity. SIM Lite will only measure the photocenter of an astrometric target, and thus determining the true astrometric orbits of such systems will be difficult. We have modified the Eclipsing Light Curve code (Orosz & Hauschildt 2000) to allow us to model the flux-weighted reflex motions of interacting binaries, in a code we call REFLUX. This code gives us sufficient flexibility to investigate nearly every configuration of interacting binary. We find that SIM Lite will be able to determine astrometric orbits for all sufficiently bright interacting binaries where the primary or secondary star dominates the luminosity. For systems where there are multiple components that comprise the spectrum in the optical bandpass accessible to SIM Lite, we find it is possible to obtain absolute masses for both components, although multi-wavelength photometry will be required to disentangle the multiple components. In all cases, SIM Lite will at least yield accurate inclinations, and provide valuable information that will allow us to begin to understand the complex evolution of mass-transferring binaries. It is critical that SIM Lite maintains a multi-wavelength capability to allow for the proper deconvolution of the astrometric orbits in multi-component systems.Comment: 12 pages, 6 figures, 6 tables. Accepted for publication in the Astrophysical Journa

Gravity and Light: Combining Gravitational Wave and Electromagnetic Observations in the 2020s

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations

A luminosity distribution for kilonovae based on short gamma-ray burst afterglows

The combined detection of a gravitational-wave signal, kilonova, and short gamma-ray burst (sGRB) from GW170817 marked a scientific breakthrough in the field of multi-messenger astronomy. But even before GW170817, there have been a number of sGRBs with possible associated kilonova detections. In this work, we re-examine these "historical" sGRB afterglows with a combination of state-of-the-art afterglow and kilonova models. This allows us to include optical/near-infrared synchrotron emission produced by the sGRB as well as ultraviolet/optical/near-infrared emission powered by the radioactive decay of $r$-process elements (i.e., the kilonova). Fitting the lightcurves, we derive the velocity and the mass distribution as well as the composition of the ejected material. The posteriors on kilonova parameters obtained from the fit were turned into distributions for the peak magnitude of the kilonova emission in different bands and the time at which this peak occurs. From the sGRB with an associated kilonova, we found that the peak magnitude in H bands falls in the range [-16.2, -13.1] ($95\%$ of confidence) and occurs within $0.8-3.6\,\rm days$ after the sGRB prompt emission. In g band instead we obtain a peak magnitude in range [-16.8, -12.3] occurring within the first $18\,\rm hr$ after the sGRB prompt. From the luminosity distributions of GW170817/AT2017gfo, kilonova candidates GRB130603B, GRB050709 and GRB060614 (with the possible inclusion of GRB150101B) and the upper limits from all the other sGRBs not associated with any kilonova detection we obtain for the first time a kilonova luminosity function in different bands.Comment: Published in MNRAS, 24 pages, 14 figure

Nuclear Physics Multimessenger Astrophysics Constraints on the Neutron Star Equation of State: Adding NICER's PSR J0740+6620 Measurement

In the past few years, new observations of neutron stars (NSs) and NS mergers have provided a wealth of data that allow one to constrain the equation of state (EOS) of nuclear matter at densities above nuclear saturation density. However, most observations were based on NSs with masses of about 1.4 M⊙, probing densities up to ∼three to four times the nuclear saturation density. Even higher densities are probed inside massive NSs such as PSR J0740+6620. Very recently, new radio observations provided an update to the mass estimate for PSR J0740+6620, and X-ray observations by the NICER and XMM telescopes constrained its radius. Based on these new measurements, we revisit our previous nuclear physics multimessenger astrophysics constraints and derive updated constraints on the EOS describing the NS interior. By combining astrophysical observations of two radio pulsars, two NICER measurements, the two gravitational-wave detections GW170817 and GW190425, detailed modeling of the kilonova AT 2017gfo, and the gamma-ray burst GRB 170817A, we are able to estimate the radius of a typical 1.4 M⊙ NS to be 11.94-0.87+0.76 km at 90% confidence. Our analysis allows us to revisit the upper bound on the maximum mass of NSs and disfavors the presence of a strong first-order phase transition from nuclear matter to exotic forms of matter, such as quark matter, inside NSs